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Isotherms of the speed of sound in helium 4 <ias as a function of pressure have been measured. 
Five of the isotherms in the liquid helium 4 temperature ran<je have been analyzed, and the results are 
discussed. The acoustically determined values of temperature are consistently higher than the values 
of temperature associated with the liquid helium 4 vapor pressure scale (T.-,h). The deviations arc 
0.006 °K at 2.323 °K, 0.007 °K at 2.807 °K, 0.008 °K at 3.211 °K, 0.010 °K at 4.212 °K and 0.012 °K at 
5.024 °K. 
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1. Introduction 

A program of absolute and secondary low tempera- 
ture thermometry has been conducted for some years 
at the National Bureau of Standards. Certain of the 
results are in the process of being reported and this 
paper presents absolute temperature determinations 
in the lower range of investigation. 

Thus far, absolute thermometry at these tempera- 
tures has been conducted by employing either a gas 
thermometer or apparatus designed for P-V isotherm 
determinations. Each of these methods is subject to 
inaccuracies which require corrections in order to ob- 
tain "the thermodynamic temperature." The correc- 
tions are caused by or related to "dead space volumes", 
deviations of the thermometric fluid from "ideal gas" 
behavior, conservation of the thermometric fluid dur- 
ing measurements, and volume changes accompany- 
ing necessary temperature changes. In some cases 
the pressure measurements, the requirements of which 
are far from being trivial, must be modified to allow 
for thermomolecular effects. To some extent, these 
corrections have been considered as existing nuisances 
but not always necessarily crucial in determining the 
'end result, a temperature. While it is not the inten- 
tion of this paper to treat the preceding corrections, the 
investigation being reported was initiated with an 
intention of avoiding them. Additionally we felt that 
a fundamentally different approach to low temperature 
thermometry was completely justified and would prove 
to be beneficial in contributing knowledge about exist- 
ing temperature determinations and scales. 

A salient feature of gas thermometry should be 
mentioned. The thermometer must be filled with the 



thermometric gas and the pressure measured at some 
assumed value of temperature. This assumption is 
critical and receives careful consideration, naturally, 
but a cascading of errors from one temperature range 
to another is possible. This type of error is also 
avoided in our approach to thermometry. 

The method of temperature determinations which 
we have pursued is based on the speed of sound in 
helium gas as a function of the gas pressure. Since 
translational dispersion of speeds in our experiments 
is negligible, one can proceed from the basic equation 

W= vElp where W is the speed of sound; 

E, the adiabatic bulk 
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modulus of the medium; and p the density, 
can be written as 
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If the equation of state for helium gas is assumed to be 
of the form 



pv = RT[\+-+^+. . . 
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the partial derivative, (t~) , can be substituted into 

eq (2) and C p /C v can be related to an ideal gas with 
the result 
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In eq (4), T is the desired temperature; 
-*> -5/3; #,, = 8.314 X 10 7 (erg/°K mol); 1 



C v /p=o 



M He = 4.0026 g/mol [2]; 2 
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anda^Kfrf+Ap^). * and C are 
the second and third virial coefficients, respectively, 
for helium 4 gas. Equation (4) suggests the form 



r 2 =r 2 u + ap+/3p 2 +. . .), 
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in which JFo is the speed of sound in an ideal gas. 
W and p are measured experimentally and treated so 
that Wo may be determined. In practice this may be 
done in two ways. W 2 and p as in eq (5) may be fitted 
by a least squares method with the intercept W%, at 
p = resulting. Another approach is to obtain 
Wo or Wq from a linear extrapolation (either graphi- 
cally or analytically). The implication is, naturally, 
that in eq (5) ap is small compared with 1, and terms 
involving higher orders of pressure negligible. Values 
of Wl obtained by the methods described afford a 
simple calculation of the temperature. Obviously, 
this treatment can produce information on helium 4 
virial coefficients, but this will be presented in a follow- 
ing paper. 

The experimentally determined quantity, speed, 
has been obtained with an apparatus which is generally 
called an acoustical interferometer [3]. This instru- 
ment employs a quartz crystal, operating at its resonant 
frequency (1 megahertz), and it is assumed that the 
crystal acts as a plane wave radiator. Pressure 
measurements are made with a mercury manometer 
which is accurate to 70 \x of mercury in the range of 
measurements. While there are a few references to 
our instrument in the literature [4, 5, 6], a detailed 
description has not been published (see footnote 3). 



2. Results 

There are listed in table 1 values of pressure, speed 
and speed squared for each isotherm. The tem- 
perature of an "isotherm" was maintained nearly 
constant by an automatic regulator which operated 
from a germanium resistor connected to a Wheatstone 
bridge. This and other germanium resistors were 
located in thermal contact with the acoustical ther- 
mometer apparatus so that allowance could be made 
for small relative temperature changes which occurred 
during measurements of an isotherm. Since four or 
five pressure-speed determinations on a given isotherm 



require a day it was not possible for all points on an 
isotherm to be taken at exactly the same temperature. 
Consequently, small deviations in temperature oc- 
curred (not exceeding 3 mdeg) and the speeds were 
corrected using the germanium thermometer indica- 
tions. The speeds in table 1 have been treated in 
this manner. 



Table 1. Isotherm data {sneed — meters/sec; pressure — 
atmospheres) 





1. 2.323 °K 




4. 4.212 °K 
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W 
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A 


0.0135 


88.991 


7.919.4 


-12.4 


0.0322 


120.241 


14.461 


-9.5 


.01541 


88.968 


7.915.3 


-1.1 


.0328 


120.182 


14.444 


-21.6 


.02005 


88.806 


7.886.5 


-7.6 


.0501 


120.030 


14.407 


4.0 


.0251 


88.528 


7.837.2 


-1.0 


.0507 


119.996 


14.391 


-1.9 


.0251 


88.527 


7.837.4 


-1.2 


.06638 


119.76 


14.342 


-2.2 


.02076 


88.740 


7.874.8 


1.6 


.0668 


119.766 


14.344 


-0.7 


.02966 


88.292 


7.795.5 


- 5.9 


.0801 


119.566 


14.287 


.6 


.02901 


88.312 


7.799.0 


-0.4 


.0997 


119.262 


14.223 


-1.6 


.03837 


87.942 


7.733.7 


2.7 


.10211 


119.232 


14.216 


-.1 










.10271 

.1204 


119.268 
118.968 


14.225 
14.145 


10.7 










2.8 




2. 2.807 ° 




.1505 


118.498 


14.042 


-0.7 










.1674 

.1986 


118.25 
117.786 


13.975 
13.874 


1.2 










3.8 


0.01601 


98.043 


9.612.5 


-8.4 


.1976 


117.76 


13.867 


- 5.9 


.0233 


97.879 


9.612.3 


3.5 


.2017 


117.707 


13.855 


-3.7 


.0301 1 


97.648 


9.535. 1 


-0.6 


(.2357) 


117.281 


13.746 




.03666 


97.445 


9.495.4 


-.7 ■ 


(.3024) 


116.139 


13.481 




.04333 


97.241 


9.455.8 


-.1 


(.2555) 
(.3805) 


116.859 
114.743 


13.646 
13.160 






3. 3.211 °K 


























5. 5.024 C K 






104.647 


10.951. 


- 5.3 




0.03048 










.04007 


104.448 


10.909.4 


3.7 


0.0644 


131.178 


17.208 


-12.5 


.050] 


104.162 


10.849.7 


- 3.0 


.0988 


130.90 


17.135 


6.2 


.0568 


103.996 


10.815.2 


- 3.3 


.1340 


130.51 


17.033 


-2.0 


.0668 


103.769 


10.768.0 


3.4 


.1676 


130.166 


16.943 


2.2 


.08366 


103.324 


10.675.8 


.2 


.1978 


129.87 


16.866 


1.2 


.0089 


102.942 


10.597. 


1.9 










.1003 


102.884 


10.585.1 


-2.6 










(.1340) 


101.995 


10.403.0 












(.1672) 


101.049 


10.210.9 












(.2036) 


99.975 


9.995.0 













1 The value of the gas constant is taken from reference [1). 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 Note added in proof: The instrument is described in Science 150, No. 3693, p. 155 
(8 Oct. 1965). 



a The columns under A contain the difference between the experimental 
W' 1 at the indicated pressures and the W' 1 of linear function determined l>\ 
the computer least squares rittin<j. 



For each isotherm, values of p and W 2 have been 
treated analytically by computer and graphically to 
obtain W%. Equation (5) demonstrates that a linear 
treatment of the data will be inadequate at suffi- 
ciently high pressures. It is not surprising, therefore, 
that in figure 1, which depicts the 3.2 °K isotherm, the 
data points markedly depart from a straight line above 
0.17 atm of pressure. In order to represent the 
isotherms linearly, we noted the departure of the 
highest point and, assuming that such departures 
are proportional to p 2 , discarded all points for which 
the deviation from linearity would exceed the experi- 
mental reproducibilities. Accordingly, the data 
above 0.1003 and 0.2017 atm for isotherms 3.2 and 4.2 
°K, respectively, were omitted in the linear computer 
analysis. For the remaining isotherms, comparable 
data selection was unnecessary since correspondingly 
high pressure data had not been acquired. 

In general, the lowest pressure points on an isotherm 
warrant lesser weightings. The inaccuracy of these 
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PRESSURE (aim) P — ► 

FIGURE 1. W 2 as a function of pressure at 3.211 °K. 

points is caused by a lack of resolution resulting from 
the high attenuation of sound energy in the thermo- 
metric gas. Consequently, the italicized data received 
one quarter weights in all computer fits. 

Column 1 of table 2 contains the W\ intercepts 
obtained from graphical treatment of the data where 
W l was plotted as a function of p. The procedure 
was straightforward in that straight lines were obvious 
representations in the graphs if the highest pressure 
points for the 3.2 and 4.2 °K isotherms were not con- 
sidered. Column 8 of table 2 contains the tempera- 
tures derived from the J^o's of column 1. 

In the computer analysis of the apparently linear 
data, each isotherm was fitted by the method of least 
squares to a polynomial in which the square of the 
speed is equated to a power series of the pressure, i.e., 
W*=Wl+Wlap + Wfip 2 + Wlyp*+. ... The first 
analysis involved the constant and linear terms, the 
second included the quadratic, and for the third the 
cubic was also added. Each of these cases yields a 
value for W\. Columns 2 and 9 of table 2 list the 



intercepts and temperatures obtained from the linear 
computer fitting. The higher order fits, with one 
exception, proved to be of little significance since the 
additional coefficients were statistically uncertain. 
The exception was the quadratic fit for the 4.2 °K 
isotherm. In this case a W% value of 14.570 with an 
uncertainty of 6 was obtained from which a tempera- 
ture of 4.209 (±0.002) °K resulted. 4 In the above 
analyses the bracketed high pressure data points for 
the 3.2 and 4.2 °K were not included. 

Equation (5) can be accurately represented by a 
linear expression (for a constant a) if the pressure is 
sufficiently low. (Our input data for the majority of 
isotherms have been selected to meet this criterion 
or, more specifically, data have been used that graphi- 
cally exhibit a linear relationship within the limits of 
the experimental error.) For this reason we believe 
that the most accurate determinations of isotherm 
temperatures are those which have resulted from a 
linear treatment. Nevertheless for completeness, it 
is necessary to extend the computer analysis to higher 
orders of pressure and include the higher pressure 
data for isotherms 3.2 and 4.2 °K. In this case the 
linear representation is inaccurate, but there exists 
no a priori approach to predict the significance of par- 
ticular higher order pressure terms since the analytical 
representation of higher order virial coefficients is 
unknown. Consequently we have approached the 
analyses empirically, by including quadratic, and 
quadratic and cubic terms for the pressure. 

When only the quadratic term is added, a statisti- 
cally meaningful coefficient for \r of eq (5) results. 
Columns 5 and 11 list the W\ and T values which result 
from extrapolations of the quadratic analyses. The 
addition of a cubic term was not helpful since the 



4 The values reported here arc not to 
mplete quadratic fitting a- reported i 



•«l with tin- value resulti 
5. 6, and 7 of table 2. 
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malysis of data 


Quadratic analysis 
of data 


Linear analysis 
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Graphical 
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Column No. 1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


13 


17.393 


17.302 


9.6 


5.4 








5.024 


5.024 


0.002 






0.012 


14.580 


14.583 


3.1 


5.0 


14.568 


6. 


7.2 


4.211 


4.212 


.001 


4.208 


0.002 


.010 


11.123 


11.117 


3.7 


3.1 


11.104 


4.4 


2.9 


3.213 


3.211 


.001 


3.207 


.001 


.008 


9.722.5 


0.717.5 


6.2 


3.1 








2.808 


2.807 


.002 






.007 


8.037.5 


8.040.5 


6.9 


4.5 








2.322 


2.323 


.002 






.006 



"The S's in columns 3 and 6 of table 2 contain the uncertainties of the coefficients. W\, in columns 2 and 5. 
''The 6's of columns 10 and 12 are the corresponding uncertainties of the temperatures in columns 9 and 11. 
' <r is the standard deviation in the least squares fitting of an isotherm. 
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statistical significance of coefficients was lost. The 
analysis was not extended for we felt that to be mean- 
ingful it should include data over the largest possible 
range of pressures. Column 13 contains the tempera- 
ture differences, r isotherm — T 58 , where T 58 is derived 
from the helium 4 vapor pressure scale [7]. These 
were obtained by calibrating germanium resistors 
directly against measured vapor pressures and com- 
paring them, in the acoustical thermometer, with 
germanium resistors which had been calibrated by 
the acoustical thermometer. In one or two cases, a 
resistor associated with the acoustical thermometer 
during an isotherm measurement had been previously 
calibrated at one or more vapor pressure temperatures. 
The temperature values of T^therm are from column 9. 
Thus our values of isotherm temperatures, column 9, 
and the indicated departures from T 5S (column 13) 
afford a comparison with the T 5 g scale as portayed in 
figure 2. Except for out data points, figure 2 is taken 
from Van Dijk [8]. In essence he has plotted values 
of absolute temperature determination, in the liquid 
helium range, relative to the Tss scale. While we are 



not prepared to present a value for the absolute accu- 
racy of our values of temperature, if one believes that 
the r 58 scale is smooth with respect to the true thermo- 
dynamic temperature scale, the plot suggests that 
perhaps the present data are the most self consistent 
of any temperature measurements that have been made 
below 5 °K. (The dispersion in temperature values 
resulting from our different methods of isotherm analy- 
ses indicates that our deduced values of temperature 
are not critically dependent upon a particular method 
of analysis.) 

The reproducibility of points on an isotherm appears 
to be within ±0.002 °K although it is obvious that most 
of the isotherms need a better definition, that is, a 
greater number of points. Nevertheless, we do feel 
that the indicated conclusions of our measurements 
are important enough to be presented now rather than 
be delayed until the complete program in this tem- 
perature range has been finished. 
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FIGURE 2. A comparison with T 58 of temperatures derived from 

thermometry and the acoustical 



gas 



isotherm measurements, 
interferometer. 
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